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Abstract 

We complete our study of the self-gravitating gas by computing the fluctua- 
tions around the saddle point solution for the three statistical ensembles (grand 
canonical, canonical and microcanonical) . Although the saddle point is the same 
for the three ensembles, the fluctuations change from one ensemble to the other. 
The zeroes of the small fluctuations determinant determine the position of the crit- 
ical points for each ensemble. This yields the domains of validity of the mean field 
approach. Only the S wave determinant exhibits critical points. Closed formulae 
for the S and P wave determinants of fluctuations are derived. The local properties 
of the self-gravitating gas in thermodynamic equilibrium are studied in detail. The 
pressure, energy density, particle density and speed of sound are computed and an- 
alyzed as functions of the position. The equation of state turns out to be locally 
p(r) = T pv(r) as for the ideal gas. Starting from the partition function of the 
self-gravitating gas, we prove in this microscopic calculation that the hydrostatic 
description yielding locally the ideal gas equation of state is exact in the N = oo 
limit. The dilute nature of the thermodynamic limit (N ~ L — > oo with N/L fixed) 
together with the long range nature of the gravitational forces play a crucial role in 
obtaining such ideal gas equation. The self-gravitating gas being inhomogeneous, 
we have PV/[NT] = f(jj) < 1 for any finite volume V. The inhomogeneous particle 
distribution in the ground state suggests a fractal distribution with Haussdorf di- 
mension D, D is slowly decreasing with increasing density, 1 < D < 3. The average 
distance between particles is computed in Monte Carlo simulations and analytically 
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in the mean field approach. A dramatic drop at the phase transition is exhibited, 
clearly illustrating the properties of the collapse. 



Contents 



1 Statistical Mechanics of the Self-Gravitating Gas 

2 Summary of the Mean Field Results 



2 

5 



3 Calculation of the Functional Determinants: the Validity of Mean Field| 6 

3~T The Grand Canonical Ense mble 8 

9 



172 The Canonical Ensemble 



3.3 The Microcanonical Ensemble 11 



4 Local Physical Magnitudes] 14 

O Particle Distribution! 14 



4.2 Average distance between particles 16 

4.3 Local energy density and gravitational potential! 16 

4.4 The pressure at a point r and the local equation of statej 21 

4.5 The speed of sound as a function of"^ 24 



5 //-dimensional generalization 



6 The Interstellar Medium 



iscussion 



8 Acknowledgements! 



28 
29 
30 
31 



A Calculation of Functional Determinants in the Spherically Symmetric] 



"Caiel 



A.l The S-wave determinant 



31 

33 

A. 2 The P-wave determinant! 34 

35 



A. 3 The D-wave and higher waves 



B Calculation of < r > and < r z > in the mean field 



35 



1 Statistical Mechanics of the Self-Gravitating Gas 

In this paper we continue our investigation of the statistical mechanics of the self-gravitating 
gas in thermodynamic equilibrium initiated in ref.0. We work in the thermodynamic 
limit which for the self-gravitating gas means the dilute limit 



A" — > oo , V — » oo , 



N 
V 1 / 3 



fixed 



where V stands for the volume of the box containing the gas. In this limit, the energy E, 
the free energy and the entropy turns to be extensive. That is, we find that they take the 
form of A" times a function of 



n 



Gm 2 N 
LT 



or £ 



EL 



Gm? N 2 



(2) 



2 



where rj and £ are intensive variables. Namely, rj and £ stay finite when N and V = L 3 
tend to infinite, rj is appropriate for the canonical ensemble and £ for the microcanonical 
ensemble. 

In paper I we derive functional integral representations for the partition function in 
each of the three statistical ensembles. These functional integrals are dominated for large 
N by a saddle point. 

When any small fluctuation around the saddle point decreases the statistical weight 
in the functional integral, the saddle point is dominating the integral and the mean field 
approach is fully valid. In that case the determinant of small fluctuations is positive. A 
negative determinant of small fluctuations indicates that some fluctuations around the 
saddle point are increasing the statistical weight in the functional integral and hence the 
saddle point does not dominate the partition function. The mean field approach cannot 
be used when the determinant of small fluctuations is negative. 

The zeroes of the small fluctuations determinant determine the position of the critical 
points for each of the three statistical ensembles. The Monte Carlo simulations for the CE 
and the MCE show that the self-gravitating gas collapses near the critical points obtained 
from mean field. 

The saddle point solution is identical for the three statistical ensembles. This is not 
the case for the fluctuations around it. The presence of constraints in the CE (on the 
number of particles) and in the MCE (on the energy and the number of particles) changes 
the functional integral over the quadratic fluctuations with respect to the GCE. 

We compute the determinant of small fluctuations around the saddle point solution 
for spherical symmetry for all three statistical ensembles. In the spherically symmetric 
case the determinant of small fluctuations is written as an infinite product over partial 
waves. The S and P wave determinants are written in closed form in terms of the saddle 
solution. The determinants for higher partial waves are computed numerically. All partial 
wave determinants are positive definite except for the S-wave. The S-wave determinant 
for each ensemble vanishes at the respective critical point. We find for spherical sym- 
metry: r]§ c = 0.797375 . . . , r$ = 2.517551 ... and r^ c = 2.03085 . . . [see fig. The 

variable i] R appropriate for a spherical symmetry is defined as i] R = G ^ N — r\ (^) ^ 3 = 
1.61199. . . 7] . 

The reason why the fluctuations are different in the three ensembles is rather simple. 
The more contraints are imposed the smaller becomes the space of fluctuations. Therefore, 
in the grand canonical ensemble (GCE) the system is more free to fluctuate and the phase 
transition takes place earlier than in the micro-canonical (MCE) and canonical ensembles 
(CE). For the same reason, the transition takes place earlier in the CE than in the MCE. 

The conclusion being that the MF correctly gives an excellent description of the ther- 
modynamic limit except near the critical points (where the small fluctuations determinant 
vanishes); the MF is valid for iV|?7 — 7] cr i t \ ^> 1. The vicinity of the critical point should 
be studied in a double scaling limit iV — > oo, rj — > r] crit . Critical exponents are re- 
ported in paper I for rj — > r]c using the mean field. These mean field results apply for 
\v ~ Vc\ "C 1 <C N\r] — r]c\ and N 3> 1. Fluctuations around mean field can be neglected 
in such a regime. 

In paper I we expressed all global physical quantities in terms of a single function 
f(r)). This function obeys a first order non-linear differential equation of first Abel's type 
pOfl . f(rj) exhibits a square-root cut at rjo, the critical point in the CE. The first Riemann 
sheet is realized both in the CE and the MCE whereas the second Riemann sheet (where 
cy < 0) is only realized in the MCE. 
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Besides the global physical quantities computed in paper I, we compute here the local 
energy density e(r), the local particle density, the local pressure and the local speed of 
sound. 

The particle distribution pv{o) proves to be inhomogeneous (except for 77 <C 1) and 
described by an universal function of 77, the geometry and the ratio f—q/R, R being the 
radial size. Both Monte Carlo simulations and the Mean Field approach show that the 
system is inhomogeneous forming a clump of size smaller than the box of volume V [see 
figs. 3 and 5 in paper I]. The particle density in the bulk behaves as pv(o) — r D ~ 3 . That 
is, the mass M(R) enclosed on a region of size R vary approximately as 

M(R) ~C R D . 

D slowly decreases from the value D = 3 for the ideal gas (77 = 0) till D = 0.98 in the 
extreme limit of the MC point taking the value 1.6 at rjc [see Table 2]. This indicates the 
presence of a fractal distribution with Haussdorf dimension D. 

Our study of the statistical mechanics of a self-gravitating system indicates that gravity 
provides a dynamical mechanism to produce fractal structures @-[[|. 

The average distance between particles monotonically decrease with 77 in the first 
sheet. The mean field and Monte Carlo are very close in the gaseous phase whereas the 
Monte Carlo simulations exhibit a spectacular drop in the average particle distance at 
the clumping transition point T (see fig.[| ). In the second sheet (only described by the 
MCE) the average particle distance increases with 77 (see fig.|3|). 

We find that the local equation of state is given by p(r) = Tpy(r). We have thus 
derived the equation of state for the self-gravitating gas. It is locally the ideal gas 
equation, but the self-gravitating gas being inhomogeneous, the pressure at the surface 
of a given volume is not equal to the temperature times the average density of particles 
in the volume. In particular, for the whole volume: PV/[NT] — f(rf) < 1 (the equality 
holds only for 77 = 0). 

Notice that we have found the local ideal gas equation of state p(f) = T pvif) for 
purely gravitational interaction between particles. Therefore, non-ideal gas equations 
of state (as often assumed and used in the literature || || |14], fL6ll ) imply the presence of 
additional non-gravitational forces. Such non-ideal equations of state appear for quantum 
gases @| . 

The energy density turns out to be an increasing function of r in the spherically 
symmetric case. The energy density is always positive on the surface whereas it is positive 
at the center for < i] R < i] R = 1.73745 . . . and negative beyond the point rj R = rj R = 
1.73745.... 

The speed of sound is computed in the mean field approach as a function of the position 
for spherical symmetry and long wavelengths. v^(r) diverges at rj R = t] R = 2.43450 ... in 
the first Riemann sheet. Just beyond this point v 2 s (r) is large and negative in the bulk 
showing the strongly unstable behaviour of the gas for such range of values of r] R . 

Moreover, we have shown the equivalence between the statistical mechanical treatment 
in the mean field approach and the hydrostatic description of the self-gravitating gas |J- 

E3- 

The success of the hydrodynamical description depends on the value of the mean free 
path (I) compared with the relevant sizes in the system. I must be <C 1. We compute the 
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ratio l/a (Knudsen number) where a is a length scale that stays fixed for N —>■ oo and 
show that l/a ~ N~ 2 . This result ensures the accuracy of the hydrodynamical description 
for large N. 

Furthermore, we have computed in this paper several physical magnitudes as functions 
of 77 and r which were not previously computed in the literature as the speed of sound, 
the energy density, the average distance between particles and we notice the presence of 
a Haussdorf dimension in the particle distribution. 

This paper is organized as follows. In section II we summarize the main results from 
on the mean field approach, in sec. Ill we compute the small fluctuations around the mean 
field theory solution. Sec. IV presents our results for space dependent quantities as the 
particle distribution, the average distance between particles, the local energy density, the 
local pressure and the local speed of sound. In section V we discuss the generalization for 
a space with any number of dimensions. In sec. VI the molecular clouds in the interstellar 
medium are discussed as a self-gravitating gas. Discussion and remarks are presented in 
section VII whereas appendix A-B contain relevant mathematical developments. 

2 Summary of the Mean Field Results 

Let us summarize here the main results of paper I in the mean field approach for spherical 
symmetry which will be used in what follows. 
The saddle point is given by 




(3) 



Here p(r) is the particle density and x(A) obeys the equation 



x"(A) + ^x'(A) + e*W=0 , x'(0)=0 



X(0) = ; . 



(4) 



x{x) is independent of r) R , and A is related to i] R through 



A x'(A) 



V 



R 



(5) 



We have in addition, 




P(l) 




(6) 



where 




(7) 



The function /a/f( ? 7 R ) obeys the Abel equation, 





For 



r] R = it follows from eq. (||) that 



Wo) = 1 . 



(9) 
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For the main physical magnitudes in the mean field approach, [that is, from the above 
saddle point and neglecting the fluctuations around it] we find: 



P V 
NT 
F-F 

NT 
S-Sp 
N 

E 

NT 



fMF(V R ) 

3[1- f MF (v R )}-V R + log f MF (v R ) 
Q[f M F(v R )- I] +V R - log f MF (v R ) 



(10) 



We find for the speed of sound squared at the surface 



s 

T 



4 + 



G fhpiri 11 ) + {l* ~ t) ford*) + \ 



(11) 



The specific heat at constant volume takes the form 

(c v )mf = 6 f MF (v R ) ~l + V R + q 7 ^f'l ; , 

2 3 Jmf{v ) - 1 

whereas for the specific heat at constant pressure we find 
(c P ) MF = 12/ mf (t7*) - - + 



(12) 



3 , 21 (//'' 2) J Midi" ■ 



6 f M F(v R ) ~V R 



(13) 



3 Calculation of the Functional Determinants: the 
Validity of Mean Field 

The mean field gives the dominant behaviour for N — > 00. We evaluate in this section 
the Gaussian functional integral of small fluctuations around the stationary points. 

As remarked in paper I, the three statistical ensembles (grand canonical, canonical 
and microcanonical) yield identical results for the saddle point. However, the small fluc- 
tuations around the saddle take different forms in each ensemble. 

Let us recall the partition functions for the three ensembles keeping quadratic fluctu- 
ations around the saddle point. In the grand canonical ensemble (see eq.(VI.29) in paper 
I), we have 



Z GC {z,T) = e 4 ^ 



{ Jo 1 d A x[\<£ \/ 2 r <t> +47T7) e<f>W]-2nr 1 log C(r?)} 



X 




1 + ° { N 



(14) 



In the canonical ensemble, we have for the coordinate partition function, 




DYdy Q e" JV s c 



Ns^[Y(.),y ] 



l + OK N 



(15) 



6 



where 



s { c\Y{-)M 



d 3 x 



Y 2 (x) r] f d 3 x d 3 y 
p(x) 2 J \x-y\ 



Y{x)Y{y)-y / d 3 xY{x) (16) 



and 



s( v R ) = 3[l-f MF (r ] H )]~r ] H + log 



3f M F(y R ) 
Air 



In the microcanonical ensemble, we have for the coordinate partition function, 

1 



N>>1 -Ns(rj) 




DY dy ^e- Ns Mcl Y (-),mwn] 
2iri 



l + O 



N 



(17) 



where, 



^ c [Y{.)m] = sP[Y(.) iyo ] - f, [ Ps (x) Y(y) 



\x - y\ 
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Fluctuations of order higher than quadratic contribute to the 1/N corrections in the 
three ensembles. 

As remarked in paper I, the functional integral for the canonical and microcanonical 



ensembles, eqs. (|15|) and (|17D , respectively, are rather close; in the last one, eq. (|17| ), there 
is an extra integration over one variable that constrains the energy. 



In the grand canonical functional integral (J1J), we have to compute the determinant 
of the operator 

L p (x, y) = 8{x - y) [- V 2 - 4tt rj e«*>] . 
In the canonical functional integral (|16"D, we find the operator 



K p (x,y) ee 



p <f>(x) 



rj 



\x - y\ 



These operators are related by 

-Vl[K p (x 1 y)e^]=L p (x,y). 

Since, 

Lo(x,y) = -5(x- y) V 2 , K (x,y) = 5(x - y) 
eq.(|TP|) can be written in abstract form as, 



(19) 



Lo K p D — L p \ 



(20) 



Here, 



D(x,y) =5{x-y) e m . 
Therefore, taking the determinant of both sides of eq.(^0|) yields, 



Dot ( 



Dot ( =^ ) r- l ' 1 ''' 1 " 17 ' 



(21) 



where we used that for any trace class operator M as D 

log DetM = Tr log M . 
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That is, 



DetD 



3 / d 3 x 4>{x) 



The functional determinants are well defined in eq.(|2~l"D by normalizing with respect to 
the vacuum values at zero density p = = 0. 

For spherically symmetric stationary points 0(r) we can expand the determinants in 
partial waves: 



log Det 



^2(21 + 1) logDet 



1=0 



where 



L> = 6(r - r>) 



# + 2i_ 1(1 + 1) 



dr 2 r dr r 2 

We evaluate these partial wave functional determinants in Appendix A. 
We expand the fluctuations in partial waves, 



Y (&) = Vi(r) Yi >m (f) 



',22) 



Lm 



where the Y^ m (f) are spherical harmonics and the c^ m arbitrary coefficients. 

The small fluctuations determinants explicitly depend on the boundary conditions 
imposed to the fluctuations Y(x) around the mean field stationary point. 

Following the arguments by Hurwitz and Katz||15||, one assumes that outside the sphere 
(no sources) the fluctuations obey the Laplace equation and therefore 



yi(.r) 



A 

1+1 



for r > 1 



Here A is some constant. Therefore, imposing continuity for Y(x) and its radial derivative 
at r = 1 yields, 

= ^ [r l+1 yi(r)] = y[{l) + (I + 1) Vl {\) = , (23) 
We impose this condition to the solutions in Appendix A. 



3.1 The Grand Canonical Ensemble 

Evaluating the Gaussian functional integral in eq.(|T4]) yields 



Z G c(z,T) 



n»i e 



-Ns GC (v R ) 



T(V R ) , n ( 1 



N 



y/Det GC {v R 

where sgc(v R ) stands for the 'effective action' at the saddle point 

1 



N 2 



(24) 



sgc(v 



K( V ) -1=2-3 fMFW 



where we used eqs.(VI.17) and (VI. 66) from paper I. 

DetGc{v R ) stands for the determinant of small fluctuations around this spherically 
symmetric saddle point and T{r] R ) for the two-loop corrections. DetGc(v R ) can be ex- 
pressed as an infinite product over the partial waves 



Det GC (^) = Det (±A = Y[[Mv R )f +1 ■ 
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where Ai(r) R ) = Det (jf V This infinite product can be properly defined using, for 
example, dimensional regularization[nj. One finds in this way that it takes a finite value. 



A (77*) = 1-±ti r f MF {r} R ) , A 1 ( V R ) = ^ r f MF ( V R ) . (25) 



A (rj R ) and Ai(rj R ) can be computed in closed form [see eqs.(|(J) and (|82"D 

We see that Aq(i] r ) > for < rj R < r] R c = 0.797375.... Therefore, the mean field 
approach breaks down for the grand canonical ensemble at rj R c . 

Notice that the determinants for all waves except the S-wave are positive definite for 

Vc > V R > 0- 

Adding the contributions from the functional determinant to the mean field results 
eq.(|T0"D in the grand canonical ensemble yields, 

3[1 - f MF (v R )} ~V R + logf MF ( V R ) + -L logDet GC (^) + O (±= 



PV fMF(v R ) + £r T 4n log Det GC (^) + O (± 



NT JMPy ' ' 6Ndr] R to ' \N 2 

= Q[fMF(v R )-n+V R -^gf MF (v R )+^(v R ^-i)^gDet G c(v R ) + 

L - 3[W^)-J] + ^i4logDet GC (^) + C»f4 



NT WM ' U 7 2 J 2iV tfy* to ^ w ' \N 2 

where we used eqs.(VI.16), (VI. 20) and (VI. 23) from paper I. 

These results correspond to include 1/N corrections in the function f(r] R ) as follows 

f(v R ) = fM F {ri R ) + lo g Det GC ( V R ) + O (± 

Eqs.([T0|)-(]i3D permit to compute the various physical quantities in terms of f(fj R ). 
Since, 

„«AlogDe tGC (, f «) W^~r - — • 

ar] u t]q C — 7] n 

jjTp, the energy and the entropy tend to minus infinity when 7] R f f] R c . This behaviour 
correctly suggests that the gas collapses for r] R j 77^. Indeed, the Monte Carlo simulations 
yield a large and negative value for -j^ in the collapsed phase (paper I). 

We want to stress that the mean field values provide excellent approximations as 
long as N\r] R c — i] R \ » 1 in the grand canonical ensemble. Namely, the mean field is 
completely reliable for large N unless r] R gets at a distance of the order V -1 from f] R c - 

3.2 The Canonical Ensemble 

We have to compute the Gaussian functional integral in eq. (|15|) 

DY dy e~ N s c m),vo) _ (27) 

J J 

where Sq [Y(.),y ] is given by eq. (|16D . The simplest way is to find a saddle point for Y(.) 
in eq. (f2~7D , that is, a solution Y(x) of the equation 

Y(x) fd 3 yY(y) , . 

V / ,! W -2/o = 0. (28) 




p(x) J \x — y\ 
9 



It is convenient to write such solution as Y(x) = yo p(x) w(x) and shift the integration 
variable in eq.(^) as follows 

Y{x)=p{x)[y w{x)+Z{x)} (29) 
where Z[x) is the new functional integration variable and w(x) is a solution of the equation 

d 3 y p(y) w(y) 



w{x) — 1] 



\x-y\ 



1 = 0. 



(30) 



(2) 

Sq\Y[.), yo] takes now the form 



sP[Y(.),y } = - [y ] 2 I d 3 xp(x) w(x) + A[Z(.)} , 



where 



A[Z(.)] = ~ / <l*.r,,(.T) Z 2 (.r) 



i] f d 3 x d 3 y 



2 J \x-y\ 



p(x) p{y) Z(x) Z{y) . (31) 



We then have, 




DYdy e- Ns cin-),yo] = j \ DZ dyQ e -f [yo? fd* X p(3) «,(*) e -NA[Z(.)] 

1 



Det(g) f<Pxp( 



(32) 



X) w(x) 



where we used eq.(pTl) and the jacobian J of the change of variables (^) has the value 

J = e!* x +&. (33) 

In the spherically symmetric case eq.(p8|) has a spherically symmetric solution w(r) 
which can be expressed in terms of the stationary point solution <p{r) as follows, 



w(r) 



2-3 V R f MF (v R ) 



2 + r 



dr 



[w(r) is related to the S-wave regular solution (|79|) 1. We can then compute the integral 
in the r. h. s. of eq.(|32|) with the result, 



d^xe^ w(r) 



2-3 V R f MF ( V R ) 
The argument of the square- root in eq.([32]) becomes then 

D C (V R ) = Det (^A f d*xp{x) w{x) = \ [3 f MF ( V R ) - l] J][A ; (^)] 2i+1 e -J*-«*) . 

(34) 

where we used eqs.(|21|), ([25|) as well as 

Ao(^) I d 3 xe^ w(r) = ± [3 f MF (r, R ) - l] 
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and we normalize to unit at r] R = 0. All factors in eq.(|34|) are positive definite except the 
first one. Hence the sign of Dc{f} R ) is defined by the sign of 3 fMF(v R ) ~ 1- 

Dc{f] R ) is thus positive for < r] R < r] R = 2.517551 . . .. That is, the mean field for 
the canonical ensemble can be applied for < r] R < rj R . 

Dc(r] R ) vanishes linearly in yr^ — r/ R at r] R = rj R . We plot the S-wave part of Dc{jj R ) 
as a function of r] R in fig. [l| 

In conclusion, the coordinate partition function e**™ in the canonical ensemble takes 
the form 



^n(v) iV > >1 



,-Ns( V R ) 



O 



For the various physical quantities, we get analogous expressions to eqs.(^), but with 
DetGc(v R ) replaced by D c {rj R ). That is, in the canonical ensemble, up to the order N^ 1 
the function f(r] R ) takes the form 



f(v R ) = fMF(v R ) + ^^og D G (r) R )+o' 1 



6N drj 1 



N 2 



(35) 



The clumping phase transition takes place when Dc{f] R ) vanishes at 1] R = rj R . Near 
such point the expansion in 1/N breaks down since the correction terms in eq. (|35|) become 
large. Mean field applies when N\r) R — r] R \ » 1. 

Since 



Vc 



12(7?* 



— oo 



(36) 



eq.(|36|) correctly suggests that PV/[NT], E/[NT] and the entropy per particle become 
large and negative for r] R f rj R . Indeed the Monte Carlo simulations yield a large and 
negative value for these three quantities in the collapsed phase (paper I). 

3.3 The Microcanonical Ensemble 

We have to compute the Gaussian functional integral in eq. (O) 




D Ydy ^-e- Ns Mcl Y (-)m,v] 
2m 



(37) 



fry\ ,w—, 

where s)J c [Y (.) , y , fj] is given by eq. (jT|) . 

As for the canonical ensemble, we start by finding a saddle point Y(.) in the Gaussian 
functional integral (|3~TD. We have, 



Y(x ) 



, 3 Y(y) + p(y) 
zsr-1] dy — ^ — ^ 1/0 = 



which has as solution 



Y(x) = p(x) 



\x - y\ 



Vo H w{x) 

VsJ Vs 



(38) 



Here, w(x) obeys eq.( 

We define a new integration variable Z{x) in eq. (|37|) as 



Y(x) = Y(x)+p(x) Z(x) 
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-0.8 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

Figure 1: The S-wave determinants A (r] R ), Dc(f] R ) and Dmc{v R ) i n the grand canoni- 
cal, canonical and microcanonical ensembles, respectively, as functions of rj R . Notice that 
the mean field approximation for each ensemble breaks down as soon as the respective 
determinant becomes negative when rj increases starting from rj = 0. 
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Eq. (p7[) takes then the form, 



DY dy ^-e- Ns Mo[n.)mM 
2ixi 




J 




DZdy ^- e - N [ A o(yo) 2 +Bov 2 +c fiyo+Alz{.)]] 



2m 



^Det(g) [AA B -Ci] 



(39) 



where 



A n = — 



d 3 x p(x) w(x) , B = 



4r? 2 77 2 



?7 2 / d 3 xp(x) w(x) , C 



d 3 xp(x) w{x) 



A[Z{.)\ is given by eq.([JT]), J by eq.(|3"3"D and we have normalized the integral to unit at 
7] = 0. 

Then, the argument of the square- root in eq. (|39"D becomes 
_ 2rf 



D M c(v 



-Det I ^ 
3 \K 



[AA Q B,-Cl] 



6 f MF (v R ) - [ ^ - v K ) fMv") + i 
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R 



HiM^)] 2 ^ 1 e~f d3xm , (40) 



J i>i 



where we used 



\ ( V R f MV R ) [4A B - Cl] = 6 f MF (v R ) " (y " V*) Jmf(v R ) + \ 



(41) 



All factors in eq.(HDf) are positive definite except the first one. The sign of Dmc(v R ) 1S 



therefore determined by the sign of the expression (41). 

Thus, Dmc{v R ) is positive in the interval < rj R < rj R and keeps positive in the 
second branch of fMF{v R ) f° r Vmc = 2.03085 . . . < rj R < rj R . At the point MC (j] R = 
rj MC = 2.03085. . . in the second Riemann sheet), the expression (|4lD becomes negative 
and the mean field approximation breaks down for the microcanonical ensemble. Dmc{v R ) 
vanishes linearly in r] R — t] MC at f] R = i] R Ic . 

In conclusion, the coordinate partition function u>(£, N) in the microcanonical ensem- 
ble takes the form 



w(£,N) 



N»l 



-Ns(r] R ) 



1 + 



^D MC (ll R ) _ 

Adding the contributions from the functional determinant to the mean field results 
([10|) yields expressions analogous to eqs.(^) but with Detcc(v R ) replaced by Dmc(v R )- 
That is, in the microcanonical ensemble the function f(f] R ) to the order iV _1 takes the 
form, 



f(v R ) = f MF (v R ) + ^gD MC (rj R ) + o(^ 



(42) 



Vmc 



QN dri R 

For 7] R I i] R IC , reaching the point MC, we find 

V R d , n / R\ rPlr&o 

T s c ' = -i2(„«-< c ) 

We see that the MF predicts that pV/[NT] grows approaching the critical point MC. 
This behaviour is confirmed by the Monte Carlo simulations. At the point MC, pV/[NT] 
increases disco ntinuously by 50% in the Monte Carlo simulations. 



+oo 
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4 Local Physical Magnitudes 



We obtain in this section physical magnitudes at at point r in the gaseous phase, that is, 
the space dependence of the particle density, the local energy density, the pressure and 
the speed of sound. We also compute the average distance between particles. 

4.1 Particle Distribution 

The particle distribution at thermal equilibrium obtained through the Monte Carlo sim- 
ulations and mean field methods is inhomogeneous both in the gaseous and condensed 
phases. In the dilute regime ij<1 the gas density is uniform, as expected. 

We plot in figs. 3-6 in paper I the density of particles from Monte Carlo simulations 
in the cube for the gaseous and for the condensed phases, respectively. 

In the mean field approximation and for the spherically symmetric case, the particle 
density is given by 

j., x A 2 e* (Ar) 

Pmf(j~) — g — ~ jj— , 0<r<l. (43) 
The mass inside a radius r is then given by 

pT \ 2 

M{r) = 4tt / r' 2 p MF {r') dr' = ^ x'(Ar) , 

Jo V 

where we used eq.(f|). For small r this gives 

M(r) ^ ^ [1 + 0(A 2 r 2 )] . (44) 

We find an uniform mass distribution near the origin. This is simply explained by the 
absence of gravitational forces at r = 0. Due to the spherically symmetry, the gravitational 
field exactly vanishes at the origin. The particles exhibit a perfect gas distribution in the 
vicinity of r = 0. Actually, eq.(|44j) is both a short distance and a weak coupling 
expression. Eq. ([44]) is valid in the dilute limit ^ R <1 for all < r < 1. 



T) 


D 


C 


0.1 


2.97 


1.0 


rj§ c = 0.797375 . . . 


2.75 


1.03 


2.0 


2.22 


1.1 


rjg = 2.517551... 


1.60 


1.07 


ri« c = 2.03085 . . . 


0.98 


1.11 



TABLE 2. The Fractal Dimension D and the proportionality coefficient C as a function 
of r) R from a fit to the mean field results according toM(r)~Cr D . 

We plot in fig. the particle distribution for 90% of the particles for several values of 
rj R . We exclude in the plots the region M(r) < 0.1 where the distribution is uniform. 
We find that these mass distributions approximately follow the power law 

M(r)~Cr D (45) 
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where, as depicted in Table 2, D slowly decreases with X(t] R ) from the value D = 3 for 
the ideal gas (rj = 0) till D = 0.98 in the extreme limit of the MC point. 

The presence of a critical region where scaling holds supports our previous work in 
the grand canonical ensemble based on field theory^, [|, ^ 

The values of the fractal dimension D of the self-gravitating gas are around D = 2 
(see Table 2). Such value can be analytically obtained assuming exact scale invariance, 
the virial theorem and the extensivity of the total energy in the limit defined by eq.(|I]) as 



follows, 18 



Let us assume that the density scales as p{r) ~ r a . Then, using the virial theorem, 
the total energy E will scale as 

E ~ f R d 3 r d 3 r' P^LP^l ~ i?5-2a 



Therefore, 



R 3 



Now, extensivity requires E/V to be independent of R, that is, a = 1 and 

p(r) ~ r- 1 , M(R) ~ R 2 . 

In addition, this implies that the gravity force in the surface of the gas is independent 
of R. 
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4.2 Average distance between particles 

We investigate here the average distance between particles < r > and the average squared 
distance < r 2 >. The study of < r > and < r 2 > as functions of 77 permit a better 
understanding of the self-gravitating gas and its phase transition in the different statistical 
ensembles. 

These average distances are defined as 



<r> = J J |r — f'| < p{r) p(r') > d 3 r d 3 r' , 
< r 2 > = / / \ f-f'\ 2 <p(f)p(f')> d 3 r d 3 r' (46) 



In the mean field approximation we have 

< p(f) p{f') >= p MF {r) p MF (f') +o(V) • (47) 

In addition, in the spherically symmetric case we use eq. (K3) for the particle density. In 
Appendix E we compute the integrals in eqs.(|46"D and we get as result, 



<r> = 2-2 / r dr 



1 r , / \ y / -1 \ ~i 2 







\ , 4>(r) - 0(1) 

T] R 



12 f 1 

<r 2 > = 2 - — / r 2 dr [0(r) - 0(1)] (46 
V Jo 



where 0(r) is given by eq.(^). 

We plot < r > and < r 2 > as functions of r] R in fig. ^. Both < r > and < r 2 > 
monotonically decrease with X(rj R ). Their values for the ideal gas are 

< r >\ n = — = 1.02857... , < r 2 >\ =-. 
177=0 35 ' lj ?=° 5 

At the critical points (C for the canonical ensemble and MC for the microcanonical en- 
semble) the average distances sharply decrease. Both < r > and < r 2 > have infinite 
slope as functions of r] R at the point C. 

We plot in fig. |4] the Monte Carlo results for < r > in a unit cube together with 
the MF results in a unit sphere. Notice that < r > sharply falls at the point T clearly 
indicating the transition to collapse. 

4.3 Local energy density and gravitational potential 

The gravitational potential at the point q is given by 

U(a) = -Gm £ _i- = £ pi- (49) 

^ — ' \q — an L L — ' r — r ; 

l<Z<iV ^ l^i^w 1 1 

where q = Lf and q] = Lri, 1 < I < N. U(q) can be easily related to the saddle point 
solution in the mean field approach. We write the sum in eq. ( |4*9"D in terms of the particle 
density p s (r) = e 9 ^ as, 

GmN f d 3 v , . 

U(q) = j- J ^ P {y) (50) 
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Figure 3: Mean value of the distance between particles (< r >) and mean value of the 
squared distance between particles (< r 2 >) as functions of rj R in the mean field approach 
from eqs.fl48|). Notice that the particles are inside a sphere of unit radius. 



Comparison of the saddle point equation (VI. 24) in paper I and yields, 

U(q) = -- [0(f) - a a ] 
and using eq.(VI.67) in paper I we recover the relation 0] 

U{$ = -- $ s (r) . 
m 



(51) 



The local density of potential energy is thus given by, 

1 NT 

e P (r) = - m p(q) U(q) = -— \<p(f) - a s \ e ( 

while the local density of kinetic energy takes the form 



3 NT 
2V 



It is easy to check that 



(r) 



U= I d 3 q e P (r) , -NT = / d 3 q e K (r) 
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Figure 4: Mean value of the distance between particles (< r >) in a unit cube from Monte 
Carlo simulations and in a unit sphere from mean field as functions of rj R . 



where U = E- 3NT/2 follows from eq.([K)]). 

In the spherically symmetric case, the local energy density takes the form 

NT \ 2 

c(r) = e K (r) + e P (r) = — [3 - V R + X (A) - X (Ar)] , (52) 

where we used eqs.(|3|), (^) and eq.(VI.55) from paper I. 
The energy density at the surface is always positive: 

e ( 1 ) = ^^( 3 -^)> ' 
V 8tti] u 

whereas the energy density at the center, 

is positive for < r] R < i] R = 1.73745. . . and negative beyond the point r] R = rj R = 
1.73745.... 

We plot in figs. [5] and ^| the energy density as a function of r for different values of r] R 
in the first and second sheets. 
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Figure 5: The local energy density, e(r) as a function of r in units of ^ for r] R = 
0.1, 1.0, 1.5 and 2.0. 
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Figure 6: The local energy density, e(r) as a function of r in units of ^y- for r] R in the 
second sheet r/ R = r/ R , 2.3, 2.15 and t)m C . 
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4.4 The pressure at a point f and the local equation of state 

The pressure p that we have obtained in sec. VII [see fig. 9 in paper I] corresponds to 
the external pressure on the gas. Let us now calculate the local pressure at a point r in 
the interior of the self-gravitating gas. We perform that computation in the mean field 
approach. 

The gravitational force is given by 

_ T N -> T N - 

F(q) = -m p(q) V q U(q) = — V r 0(r) = — V r e^ , (53) 

where we used eq. (pTj) , the expression for the particle density in the volume L 3 , P v(q) = 
n e 4>(r) anc [ q = Z,<r. 

Since the density of force eq.(|53|) is the gradient of the local pressure, we find 



TN -> 

VrP(r) = -jj- V r e^ = T V rPv {r) (54) 

i. e. 

p(r) = T p v (r) (55) 

That is, we have shown that the equation of state for the self-gravitating gas is locally 
the ideal gas equation in the mean field approximation. Notice that contrary to ideal 
gases, the density here is never uniform in thermal equilibrium. Therefore, in general the 
pressure at the surface of a given volume is not equal to the temperature times the average 
density of particles in the volume. In particular, for the whole volume, PV/[NT] < 1 
(except for r\ = 0). 

The local pressure in the spherically symmetric case can be written in a more explicit 
way using eqs.@ and fl55l): 

NT - 3t/« [ } 

For r = 1, eqs. (0) and ( [TDD show that p = p(l) coincides with the external pressure. 
The local density and the local pressure monotonically decreases with r. 
The particle density at the origin follows from eqs. (§3) and (P|): 



4 7T 7] 1 



P{r = 0) = 7— r (5T) 



This particle density at the origin grows when moving from rj = to MC as shown in 
fig. |7[ In particular, for i] R = we have 

p(r = 0)\ vR=0 = A (58) 

where we used eqs. (^), (Qj and (^). Notice that p(r) is r-independent for r] R = 0. 

The particle density at the surface is proportional to fMF{j] R ) [ see ec L-(H)] an d plotted 
in fig. 9 of paper I. We see that it decreases when moving from 77 = to 77 = rjMin = 
2.20731 ... in the second sheet. The migration of particles towards the center as rj varies 
is manifestly responsible for these variations in the density. 

The pressure (and density) contrast is given by 

PM = = e -x(A) 
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We plot in fig. § x(^(v R )) as a function of rj R . For r] R = rj R and r] R = 1] R IC we recover 
the known values p(0)/p(l) = 32.125 . . . and p(0)/p(l) = 708.63 . . ., respectively j]!!] , |T3 



Notice that pV/[iVT] < 1 (see fig. 9 in paper I) for rj > although the equation of 
state is locally the one of an ideal gas as we have showed. The inhomogeneous particle 
distribution in the self-gravitating gas is responsible of such inequality. 

Local equations of state other than the ideal gas are often assumed in the context of 
self-gravitating fluids [5], [TT], |TJ[]. Our result imply that forces other than gravitational 
are necessary to obtain a non-ideal local equation of state in thermal equilibrium. 

We have thus shown the equivalence for the self-gravitating gas between the statistical 
mechanical treatment in the mean field approach with the hydrostatic description p|-|16||. 

4.5 The speed of sound as a function of r 

For very short wavelengths \ s <C L, the sound waves just feel the local equation of state 
(|55|) and the speed of sound will be that of an ideal gas. For long wavelengths (of the 
order L), the situation changes. The calculation in eq.([Ll]) corresponds to the speed of 
sound for an external wave arriving on the sphere in the long wavelength limit. Let us 
now make the analogous calculation for a wave reaching the point q inside the gas. 
Our starting point is again eq.(III.29) in paper I, 

where cp and cy are the specific heats of the whole system at constant (external) pressure 
and volume, respectively, and p(q) is the local pressure at the point q. 
We find for the spherically symmetrical case in MF 

^1 = C -L [ 6 f( n R ) + Ar y'(Xr)] e x{Xr) (60) 

T c v 9 V R [3/(77*) -1] 1 J{V ) + Ar *^ r )\ ' ^ 



where we used eqs.(|56|), (|59]) and 



dr] 



dV) T 3V 



[X is a function of r] R as defined by eq.(|])]. 

At the surface, (r = 1), Vg(r) reduces to eq. (pTT|) after using eqs.(|5|), (0), (|12]) and flllf ). 



For i] R < i] R = 2.43450 ... in the first sheet, v 2 s {r) is positive and decreases with r as 
shown in fig. |9]. 

At r] R = r/ R , Vg(r) diverges for all < r < 1 due to the factor cp in eq. (|60|) [cfr. 
eq. (|13D1. The derivative of p with respect to V is proportional to 6 f(r] R ) + A r x'(A r) as 
we see in eq.(^). At r = 1 this factor becomes 6/(77 ) — V R [ see ec L-(il)] which exactly 
vanishes at i] R = rj R canceling the singularity that cp possesses at such point [see eq.([L3|)]. 
Thus, Vg(l) is regular at r] R = rj R . 

Vg(r) becomes large and positive below and near t] R and large and negative 
above and near 77^ = 2.43450 ... as depicted in fig. |10[ This singular behaviour witness 
the appearance of strong instabilities at 77 = 77^. For 77* > 77^, v s (r) becomes imagi- 
nary indicating the exponential growth of disturbances in the gas. This phenomenon is 
especially dramatic in the denser regions (the core). 
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Figure 9: v 2 (r)/v 2 (l) as a function of r for r] R = 2.0, 1.5, 1.0 and 0.1. That is, values of 
i] R smaller than r} R = 2.43450 . . .. v 2 {r)/v 2 {l) is here always positive and decreases with 
r. 
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Figure 10: v 2 s (r) / v 2 s (l) as a function of r for values of rj R around 7] R = 2.43450 

Positive values of v 2 (r)/v 2 (l) correspond to rj R < r] R and negative values of v 2 s (r) / v 2 s (l) 
correspond to i] R > i] R . We see that the speed of sound squared tends to +oo in the bulk 
(r < 1) for i] R | rj R while it tends to — oo for i] R j 77^. 
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Figure 11: v 2 (r)/v 2 (l) as a function of r for r] R in the second sheet: r] R = rj^ c = 
2.03085 . . . , 2.15, 2.3 and rjg = 2.517551 . . .. v 2 (r)/v 2 (l) is here strongly negative in the 
core of the sphere. Notice that v 2 (l) < for r] R = ?7^- c = 2.03085 . . . while v 2 (l) > for 
r] R = 2.15, 2.3 and r^. 
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For r) R beyond rj R and before rj R = 2.14675 ... in the second sheet v 2 {r) stays negative 
around the core while it becomes positive in the external regions as depicted in fig. |TT| . 
For example, v 2 (r) is positive at 7] R = rj R for r > 0.4526 .... 

For r] R in the second sheet beyond r) R and before r] R IC , v 2 (r) is positive in the core 
and negative outside. 



5 z^-dimensional generalization 



The self-gravitating gas can be studied in //-dimensional space where the Hamiltonian 
takes the form 



A? 



H 



N 



1=1 



^ 2m ^ 

l<l<j<N 



_ ->\v—2i 



and 



N 



H 



N 



y — — Gm 2 y to t 

^ 2m ^ 



1 



I Si - SiU' 



for v ^ 2 



at v — 2 . 



(61) 



(62) 



1=1 l<l<j<N 

The partition function in the microcanonical, canonical and grand canonical ensembles 
takes forms analogous to eqs.(II.2), (III. 1 ) and (VI. 7) in paper I, respectively. 
We now find for the microcanonical ensemble, 



S(E,N)= log 



jyuN-2 m 3zWV/2-2 jjv(2-v/2) N+v-2 QvN/2 



-1 



NW^f) (2tc) uN/2 
where the coordinate partition function takes now the form 



+ log w(£,N) . 



w(£,N) = 
with 



1 N 



i=i 



n 



r N ) 



-i vN/2-l 



r N ) 



EL 



u-2 



and 



Gm 2 N 2 

1 E 



(63) 



N 



u[n, ...,r N 

\ 

l<l<j<N 

In the canonical ensemble we obtain now, 

uN 



u-2 



Z C (N,T) 



1 fmTL 2 



N\ 



2tt 



1 N 



IP 

1=1 



Ti e 



r\ u(f 1 ,...,fjvr) 



where the variable r\ takes the form 



Gm 2 N 



As we can see from eqs. 
the exponent of L. 



V = • (64) 

the only change going off three dimensional space is in 



In //-dimensional space the thermodynamic limit is defined as V, N — > oo keeping 77 
and £ fixed. That is, N/L u ~ 2 = N/V 1 ^ 2 ^ is kept fixed. The volume density of particles 
N/V vanishes as V~ 2 I V for V, N — > 00 and v > 2. It is a dilute limit for v > 2. 

When v < 2, one has to assume that the temperature tends to infinity in the thermo- 
dynamic limit in order to keep rj and £ fixed as V, N — > 00. 
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6 The Interstellar Medium 



The interstellar medium (ISM) is a gas essentially formed by atomic (HI) and molecular 
(H 2 ) hydrogen, distributed in cold (T ~ 5 — 50K) clouds, in a very inhomogeneous and 
fragmented structure. These clouds are confined in the galactic plane and in particular 
along the spiral arms. They are distributed in a hierarchy of structures, of observed 
masses from 10~ 2 M to 1O 6 M . The morphology and kinematics of these structures 
are traced by radio astronomical observations of the HI hyper fine line at the wavelength 
of 21cm, and of the rotational lines of the CO molecule (the fundamental line being at 
2.6mm in wavelength), and many other less abundant molecules. Structures have been 
measured directly in emission from O.Olpc to lOOpc, and there is some evidence in VLBI 
(very long based interferometry) HI absorption of structures as low as 10~ 4 pc = 20 AU 
(3 10 14 cm). The mean density of structures is roughly inversely proportional to their 
sizes, and vary between 10 and 10 5 atoms jew? (significantly above the mean density of 
the ISM which is about 0.1 atoms/cm? or 1.6 10~ 25 g/cm 3 ). Observations of the ISM 
revealed remarkable relations between the mass, the radius and velocity dispersion of the 



various regions, as first noticed by Larson [Ell], and since then confirmed by many other 



independent observations (see for example ref.|22||). From a compilation of well established 
samples of data for many different types of molecular clouds of maximum linear dimension 
(size) R, total mass M and internal velocity dispersion Av in each region: 

M(R) ~ R d « , Av ~ R q , (65) 

over a large range of cloud sizes, with 10~ 4 — 10~ 2 pc < R < 100 pc, 

1-4 <d H <2, 0.3 < q < 0.6 . (66) 

These scaling relations indicate a hierarchical structure for the molecular clouds which is 
independent of the scale over the above cited range; above 100 pc in size, corresponding 
to giant molecular clouds, larger structures will be destroyed by galactic shear. 

These relations appear to be universal, the exponents dn, q are almost constant 
over all scales of the Galaxy, and whatever be the observed molecule or element. These 
properties of interstellar cold gas are supported first at all from observations (and for 
many different tracers of cloud structures: dark globules using 13 CO, since the more 
abundant isotopic species 12 CO is highly optically thick, dark cloud cores using HCN or 
CS as density tracers, giant molecular clouds using 12 CO, HI to trace more diffuse gas, 
and even cold dust emission in the far- infrared) . Nearby molecular clouds are observed 
to be fragmented and self-similar in projection over a range of scales and densities of at 
least 10 4 , and perhaps up to 10 6 . 

The physical origin as well as the interpretation of the scaling relations ( |65D have been 
the subject of many proposals. It is not our aim here to account for all the proposed 
models of the ISM and we refer the reader to refs. [EJ| for a review. 



The physics of the ISM is complex, especially when we consider the violent pertur- 
bations brought by star formation. Energy is then poured into the ISM either mechani- 
cally through supernovae explosions, stellar winds, bipolar gas flows, etc.. or radiatively 
through star light, heating or ionizing the medium, directly or through heated dust. Rel- 
ative velocities between the various fragments of the ISM exceed their internal thermal 
speeds, shock fronts develop and are highly dissipative; radiative cooling is very efficient, 
so that globally the ISM might be considered isothermal on large-scales. Whatever the 
diversity of the processes, the universality of the scaling relations suggests a common 
mechanism underlying the physics. 
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We proposed that self-gravity is the main force at the origin of the structures, that 
can be perturbed locally by heating sources^ |J. Observations are compatible with 
virialised structures at all scales. Moreover, it has been suggested that the molecular 
clouds ensemble is in isothermal equilibrium with the cosmic background radiation at 
T ~ 3K in the outer parts of galaxies, devoid of any star and heating sources J7|. This 
colder isothermal medium might represent the ideal frame to understand the role of self- 
gravity in shaping the hierarchical structures. 

In order to compare the properties of the self-gravitating gas with the ISM it is con- 
venient to express m, T and L in rj in appropriate units. We find from eq.(Q) 

V = 0.52193 , 

where m is in multiples of the hydrogen atom mass, T in Kelvin, L in parsecs and A4 Q 
is the mass of the cloud in units of solar masses. Notice that L is many times (~ 10) the 
size of the cloud. 

The observed parameters of the ISM clouds|22j yield an rj around ~ 2.0 for clouds 



not too large: M. Q < 1000. Such r\ is in the range where the self-gravitating gas exhibits 
scaling behaviour. 

We conclude that the self-gravitating gas in thermal equilibrium well describe the 
observed fractal structures and the scaling relations in the ISM clouds [see, for example 
fig. and table 2]. Hence, self-gravity accounts for the structures in the ISM. 



7 Discussion 

We have presented in paper I and here a set of new results for the self-gravitating thermal 
gas obtained by Monte Carlo and analytic methods. They provide a complete picture for 
the thermal self-gravitating gas. 

Starting from the partition function of the self-gravitating gas, we have proved from a 
microscopic calculation that the local equation of state p(r) —Tpv(f) and the hydrostatic 
description are exact. Indeed, the dilute nature of the thermodynamic limit (N ~ L — > oo 
with N/L fixed) together with the long range nature of the gravitational forces play a 
crucial role in the obtention of such ideal gas equation of state. 

More generally, one can investigate whether a hydro dynamical description will apply 
for a self-gravitating gas. One has then to estimate the mean free path (I) for the particles 
and compare it with the relevant scales a in the system JL9|. We have, 



1 L 3 

I ^_ (67) 

p v a t N cr t 

where py = -p- p is the volume density of particles and a t the total transport cross section. 

Due to the long range nature of the gravitational force, o t diverges logarithmically for 
small angles. On a finite volume the impact parameter is bounded by L and the smaller 
scattering angle is of the order of 

Aq G m 2 
q LT 

since q = mv ~ Vm T and Aq ~ ^£ ~ T7F' 



30 



We then have for the transport cross section |jl9fl, 



(Gm) 2 , LT 
log 



u — v 



714 



Gm 2 



77 



1 ^ 
log — 

77 



(6? 



where we used that v ~ 
large impact parameters 
From eqs. 



and 



/ 



— and \v — v'\ ~ Aq/m. As we see, the collisions with very 
L) dominate the cross-section, 
we find for the mean free path: 

2 1 T. «2 

(69) 



L 

iV 



Gm 2 



TL 



rj 



log 



T L 

Gm 2 



N 3 log 



iV 



where we have here replaced py by -p-. A more accurate estimate introduces the factor 
p = in the denominator of I. This factor for spherical symmetry can vary up to two 
orders of magnitude [see fig. [7[] but it does not change essentially the estimate 



We see from eq. (|69|) that in the thermodynamic limit I becomes extremely small com- 
pared with any length a = O(N ) that stays fixed for A" — > 00. We find from eq.( 



N 2 log £ 



In conclusion, the smallness of the ratio l/a (Knudsen number) guarantees that the hy- 
drodynamical description for a self-gravitating fluid becomes exact in the N, L — » 00 
limit for all scales ranging from the order L° till the order L. 

It must me noticed that the time between two collisions t co i — l/v^l is different 
both from the relaxation time and from the crossing time used in the literature. In 
particular, it is well known that! 



MM 



t 



crossing 



I log AT. 



^relaxation 

This formula does not concern the time t co \ between two successive collisions. The time 
t co i is indeed very short due to the small angle behaviour of the gravitational cross section. 
For constant cross sections one finds a very different result for t co i [see ref. [16[] . 



In this paper and in its companion paper |IJ we thoroughly investigate the physics of 
the self-gravitating gas in thermal equilibrium. It is natural to study now the hydrody- 
namics of the self-gravitating fluid using p(f) = T p{f) as local equation of state. A first 



work on this direction is ref. [23]. 
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A Calculation of Functional Determinants in the Spher- 
ically Symmetric Case 

We compute here the determinant of the one-dimensional differential operator: 

ar z r dr r z 
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where <p(r) is the stationary point given by eqs.(|3])-(fD. 
It is convenient to change the variable r to 

x = logr , < r < 1 , — oo < x < , 

and perform a similarity transformation by ^/r = e x l 2 on Di(l). That is, we define 

d 2 



D(l) = e~ x/2 D 1 {l)e 



x/2 



dx 2 



k +v(x) 



where, 



v{x) = -4nr] R e ^ r=e ^ +2x , k = i(l + ^\ . 
D(l) has the form of a standard Schrodinger operator. Notice that, 
lim v(x) = and v{0) = -3r] R f M F{v R ) < • 



X— > — oo 



We thus have an 'attractive' potential v(x). The appearance of a 'bound state' (that is, a 
negative eigenvalue k 2 < 0) corresponds here to an instability in the self-gravitating gas. 
We compute now the determinant 



A, = A(k) 



Det 



. d 2 
dx 2 



k 2 + v{x) 



Det [-£ - k 2 } 



where we normalize at v — 0, as usual. 

The logarithmic derivative with respect to k 2 can be expressed in terms of the inverse 
of the operator D{1) 



d log A (k) 

dk 2 



dx 



1 



1 



x 



■£-k 2 + v(x) 



\x) — \x\ 



dx 2 



\X) 



(70) 



In order to express the Green function (inverse of the operator D(l)) we introduce the 
Jost and regular solutions of the Schrodinger-like equation 



d 2 



dx 2 



— k + v(x) 



f(k,x) = 0. 



The solution of eq. (ffTD with the asymptotic behaviour 

f(k,x 



x—^_ oo ^—ikx 



(71) 



(72) 



is called the Jost solution, while the solutions ip(k,x) of eq.([7lD obeying the boundary 
conditions 

ip(k,0) = a , (p(k,0) = /3 (73) 

where a and (3 are fc-independent arbitrary parameters, are called regular solutions. The 
values of a and (3 must be selected on physical grounds as discussed in sec. III. A(/c) 
explicitly depends on them as we shall see below. 
For large and negative x we have 



(p{k,x] 



x— >— oo 



F(k) e ikx + G(k) e 



-ikx 
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where the coefficients F(k) and G(k) depend on v (.). F(k) is called the Jost function. 
The small fluctuations defined by eq.(^) are related to the regular solutions by 

yi(r) = —j= <f(k = i ( I + - ) , x = logr) (74) 



where yi(r) is the radial part of Y(x) in the Zth. wave [see eq.(p2|)]. 
The Green function of D(l) takes then the form 

1 i 
( x l ~1p TV, — rT = ttt f(k, x>) f(k, x<) , (75) 

where x> (x<) stands for the larger (the smaller) between x and x'. 

Inserting eq.(|75]) into eq . (|70|) , the integration over x can be performed with the help 
of wronskian identities with the result 

_ ^F(k) 
A[k) ~ ak-i(3> 

where we used that A(oo) = 1 as normalization. 

The Jost function F(k) can be expressed in terms of the Jost solution at the origin 
just computing at the origin the wronskian of the regular and Jost solutions, 

F(k) = ^ [(3f(k, 0) - af(k, 0) 
Therefore, once the Jost solution is known, the calculation of the determinant is immediate 

= W(M)-./M (76) 

v ' (3 + iak v ' 

Notice that A(k) is an homogeneous function of a and j3. 

We impose the physical boundary conditions ( p3|) to the regular solution (|73|). This 
gives 

(2l + l)a + 2f3 = (77) 

and 

A, l /(M) + 2[Ti^ 0) (78) 



A.l The S-wave determinant 

For / = 0, k = i/2, eg. (|7l|) can be solved in terms of the stationary point solution 
0(r = e x ). That is, 

1 d(fi(r = e x 



f(z/2,x) = e*/ 2 



1 + 



dx 



(79) 



obeys both eq.(|7T|) and the boundary conditions (|72|). Eq.(^T|) for the function ( [79|) can 
be checked just taking the derivative of the stationary point equation (VI. 43) in paper I 
with respect to x — logr. 

We can now compute f(i/2, 0) and /(i/2, 0) using eqs.(VI.43), (VI. 46) in paper I, and 
(ED with the result, 



/(i/2,0) = l- 



n 



R 



/(</2,0) 
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Figure 12: The partial wave determinants for I = 1, 2, 3 and 4 as a function of r/ fl . Notice 
that all these determinants are positive definite for t]q > r] R > 0. 



This yields for the S-wave determinant (W> 



A = A(k = z/2) 



a -2/3 



a 



for arbitrary values of a and 

For the boundary conditions ( [77D we find 



Ao(V 



1- 2^/^(1^) 



*0) 



A. 2 The P-wave determinant 



/(3i/2,a0 = -— e" 



(81) 



For 1 = 1, k = 3i/2, eg. (|71~D can also be solved in terms of the stationary point solution 
<f>(r = e x ). That is, 

_3_ e _ x/2 d(j){r = e x ) 
X 2 dx 

obeys both eq. (|7lD and the boundary conditions (|72|) . Eq.([H]) for the function (jSlD can 
be checked just taking the derivative with respect to r of the stationary point equation 
(VI. 43) in paper I. 

We obtain for /(3i/2,0) and /(3i/2,0), 



/(3i/2,0) 



3?7 J 



/(3z/2,0) = ^ 
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where we used again eqs.(VI.43), (VI. 46) in paper I, and (|). 

We finally obtain for the P-wave determinant for arbitrary values of a and (3, 



Ai = A(k = 3i/2) 



3 rj 



R 



X 2 (3a-2(3) 
For the boundary conditions (|77|) we find 



{3a[2f MF (r ] R )-l]-2f3} 



Ai(r/ 



3rj R 



fMF(V R ) ■ 



(82) 



A. 3 The D-wave and higher waves 

We plot in fig. |l2| the determinants for I = 2,3,4 and 5 for the boundary conditions 
([T?|) We see that these partial waves determinants are positive for all positive values of 

T] R < f] R . 

We can evaluate the Jost solutions asymptotically for large k (large I). Using the 
standard Riccati transformations yields 



f{k,x] 
where k — i (l + ^) and 



1 + I^L _ ^) + 2^) 
2/ + 1 2(2/ + 1) 2 



1 



(2i + 1; 



J(x) = / dyv(y) = \e x X '(^) + x(Xe x ) 

J—oo 

Therefore, for the boundary conditions (TTTf) we find from eq.([T 



a / ^ 1 J (0) / 2 (0) - / 1 

A,(r? R ) = H — — h C 

iyi ' 21 + 1 2(2/ + 1) 2 V (2Z + 1 



where, 



J(0) = log 



3 77 



R 



-fMF(V R ) 



B Calculation of < r > and < r 2 > in the mean field 



In the mean field approach we have from eqs.fllq) and (|47D 

< r >= / r 2 dr r' 2 dr' j cfO(f) / dO(r') p(r) p(r') \J r 2 + r' 2 — 2rr' cos(r, r') 



< r 2 >= / r 2 dr r' 2 dr' / dfi(f) / dVt(r) p(r) p(r') [r 2 + r /2 — 2rr' cos(r,f')] . 
Jo Jo J 

Integrating over the angles yield, 



16 7T 2 f 1 , . , f T /2 J . . ,, 
< r >= / rc/r p(r) < r dr p{r ) 

3 Jo lio 



r' 2 + 3r 2 



+ r r dr' p(r') 



r 2 + 3 r' 2 



< r 2 >= 87r / r 4 dr p(r) 
Jo 
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We have then to compute integrals of the type J r n dr e^ r > for n — 1, 2, 3 and 4. We 
find using eqs.(VI.43), (VI.45), (VI.46) and (VI.54) in paper I, 



4 7T T) R 


ire* 
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dr' = 
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_ r 2 d(f>{r) 
dr 

— [r (j){r)\ 



dr 

_ r 4#W +2r 3^ r ) _ g / r ' 2 0( r ') rfr ' ( 83 ) 
dr 



Collecting all terms yields after calculation eq. (|48|) . 
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